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Abstract. The package Binomials contains implementations of specialized algorithms for 
binomial ideals, including primary decomposition into binomial ideals. The current imple- 
mentation works in characteristic zero. Primary decomposition is restricted to binomial 
ideals with trivial coefficients to avoid computations over the algebraic numbers. The basic 
ideas of the algorithms go back to Eisenbud and Sturmfels' seminal paper on the subject. 
Two recent improvements of the algorithms are discussed and examples are presented. 



1. Binomial ideals 

Let S = k[xi, . . . , Xn] denote the standard polynomial ring over a field k. A binomial ideal 
/ C S is an ideal generated by binomials — Ax^, where u,v (z N" are exponent vectors and 
A G k is a coefficient. Monomials are also considered binomials. Assumptions on k will be 
forced upon us when computing primary decompositions. The ideal {x^ — 1) has no primary 
decomposition into binomial ideals when k does not contain a third root of unity. Interest in 
binomial ideals is due to the frequency with which they arise in applications. To name one, 
in algebraic statistics one is interested in primary decompositions of conditional independence 
ideals whose components describe various combinatorial ways in which a set of conditional 
independence statements can be realized [U [5]. Because the minimal primes of binomial 
ideals are toric ideals [2] , binomial conditional independence models are unions of exponential 
families. In particular they are unirational. Knowledge of a primary decomposition also gives 
a piecewise parameterization of such models. 

The new Macaulay 2 package Binomials offers specialized implementations of primary 
decomposition, radical computations and minimal and associated primes. The starting point 
for this implementation was Section 9 in Eisenbud and Sturmfels' foundational paper [2], but 
various improvements have been discovered and implemented. Binomials is the fastest and 
often only way to compute large primary decompositions of binomial ideals. 

Example 1. 



il 
i2 
i3 



needsPackage "Binomials" 
R = QQ[x,y] 

I = ideal (x"2-x*y, x*y-y~2) 



14 : binomialPrimaryDecomposition I 
[...] 

2 

o4 = {ideal (x - y) , ideal (x, y )} 

A binomial primary decomposition starts with a cellular decomposition. Recall that a 
binomial ideal / C S is cellular if in S/I every monomial is either regular (i.e. a nonzerodivisor) 
or nilpotent. The implemented algorithm to compute a cellular decomposition is discussed 
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in [6l[8]. Since cellular decomposition is independent of k, it can serve as a first approximation 
of primary decomposition over any field. 

In this paper we focus on decomposing a cellular binomial ideal further. To this end, assume 
that I is J-cellular for some J C [n], that is, the variables with indices in J are regular, while 
the variables with indices in J := [n] \ J are nilpotent. 

2. Computing associated primes 

If k is algebraically closed, then the associated primes of a binomial ideal are guaranteed 
to be binomial. Since computer algebra system usually don't implement algebraically closed 
fields, the input binomial ideals are restricted to be generated by unital binomials — x'" . 
In this case the binomial primary decomposition together with the associated primes exist 
over a cyclotomic extension of Q [6] . If necessary. Binomials will construct this extension and 
return its result over a different ring. 

Example 2. 

il : R = QQ[x] 
±2:1= ideal(x-3-l) 
i3 : BPD I 
[...] 

o3 = {ideal (x - 1), ideal (x - ww ) , ideal(x + ww + 1)} 

3 3 

In the following discussion we will assume k to be algebraically closed and of characteristic 
zero. Let / be J-cellular for J C [n], and denote mj := {xi : i ^ J). The associated primes 
of / are of the form Ip^j + mj, where Ipj := {x'^ — p{u — v)x'" : u — v G L) is a lattice ideal 
in the J-variables, and /? : L — )> k* is a group homomorphism from a sublattice L C Z'^. The 
pair {p,L) is called a partial character in [U |2]. Here we will simply speak of a character. 
An extension of p is a character r : L' — )• k such that L Q L' and p and r agree on L. A 
character is saturated if its domain is saturated lattice, that is L is not contained with finite 
index in any other sublattice of Z"^. An extension to a saturated character is a saturation. 
Denote k[J] := k.[xi : i £ J], and k[J] := k[xi : i ^ J]. Associated primes of cellular binomial 
ideals come in groups. The following theorem states that they are to be found among the 
associated primes of lattice ideals in k[J]. 

Theorem 3 (O Theorem 8.1]). Let I S be a J-cellular binomial ideal. Let I^^j + mj be 
an associated prime of I, then there exists a monomial m S k[J], and a character r on Z"^ 
whose saturation is a, such that 

{I ■.m)n k[J] = Jr. 

It can be seen that the converse also holds. Every associated prime of any occurring lattice 
ideal is associated to I [7]. Theorem [3] shows that a sub-problem in the computation of 
associated primes is to determine the set of lattice ideals of the form (/ : m) n k[J]. 

Definition 4 (H)- A lattice L C Z"' is potentially associated to / if there exists a witness 
monomial m G k[J] such that (/ : m) PI k[J] = Ip^j for some character p : L — )■ k*. 

Properly defining the set of associated lattices, which is contained in the set of potentially 
associated lattices, requires care and is one of the topics of [7]. For computational purposes 
the few lattices that are potentially associated, but not associated, play a minor role. They 
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will eventually yield redundant primary decompositions, a problem that has to be handled 
in any case, since cellular decomposition introduces redundancy on a large scale. The lattice 
ideals (I : m) n k[J] are partially ordered by inclusion, and so are their lattices. 

Definition 5. A potentially associated lattice is called embedded, if it properly contains the 
lattice of /nk[J]. 

A first algorithm to find potentially associated lattices would examine all ideals (/ : m) 
where m is a nonzero monomial in k[J]/(I n Ik[J]). By cellularity of / there are only finitely 
many such monomials and this search will terminate. The associated primes algorithm in 
Binomials instead uses a random search. The set of monomials to be examined can be very 
large compared to relatively few potentially associated lattices. The design goal in Algorithm [1] 
is to compute as few colon ideals (/ : m) as possible. If a monomial m divides a monomial 
n, then (/ : m) C (/ : n) and containment also holds for the potentially associated lattices. 
Due to this fact we can exclude large posets of monomials if we find two monomials with the 
same potentially associated lattice. 

Algorithm 1. 

Input: A J -cellular binomial ideal I. 

Output: The potentially associated lattices of I 

(1) Compute the lattice ideal I PI k[J]. 

(2) Initialize a list of known potentially associated lattices and witnesses containing only 
the pair (In k[J], 1). 

(3) Initialize a todo-list with all monomials in a h-basis o/k[J]/(/ fl k[J]). 

(4) Iterate the following until the todo-list is empty 

• Choose and remove a random monomial m from the todo-list. Compute the lattice 
ideal (/ : m) fl k[J] and check if its lattice is already on the list of potentially 
associated lattices. 

— If yes, then add m as a new witness for that lattice, remove from the todo- 
list every monomial between existing witnesses and m. 

— If no, then add (/nk[J],m) to the list of potentially associated lattices. 

To save space and time, the implementation in Binomials does not save all the witness 
monomials. If m, n are both witnesses for the same potentially associated lattice and m|n, 
then only m needs to saved. 

Given the set of potentially associated lattices, determining the associated primes is easy. 
It consists of saturating characters and will not be discussed here. The necessary cyclotomic 
extensions are handled in a separate package Cyclotomic, published together with Binomials. 

3. Computing minimal primary components 

Let / C 5* be a J-cellular binomial ideal, and P = Ip^j + mj one of its associated primes. 
Eisenbud and Sturmfels show that any primary component of / over P contains Ipj. In fact, 
I-\-Ip^j has P as its unique minimal prime and a primary component is computed by removing 
all embedded primary components from I -\-Ip^j. This is the content of [21 Theorem 7.1]. Let 
Hull(/) denote the intersection of the minimal primary components of a binomial ideal /. If 
/ is cellular, then Hull(/) is binomial. Computing Hull of a binomial ideal is a cumbersome 
procedure. One way, described in [2], is to successively identify binomials h such that (I : h) 
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is a binomial ideal strictly containing /. This approach is slow. Here we will use a similar 
strategy like in Algorithm [TJ Denote McmbC-^) the monomial ideal generated by all witnesses 
of embedded lattices of /. Then [H Theorem 3.2] implies the following simplification. 

Proposition 6. // 1 is J -cellular and has exactly one minimal prime, then 

Hull(/) = / + Me^b(/) 

In particular Hull(/) is binomial. 

To compute the minimal primary component of / over P = /p^j+mj one computes Hull(/+ 
Ip,j) [21 [8] • The monomial ideal Memb(-^ + -^p,j) is determined essentially by Algorithm [TJ It 
is in fact simpler, since only minimal generators of Meinb(-^) need to be computed. In most 
cases only a small fraction of the standard monomials needs to be examined. 

Example 7. This example demonstrates how a component of high multiplicity leads to many 
monomials to be examined. 

11 : R = QQ[a,b] 

12 : I = ideal (a'lOOOO * (b-1)) 

13 : BPD I 

10000 

o3 = {ideaKb - 1), ideal(a )} 

In this case the poset of nilpotent monomials is totally ordered and a typical run of Algorithm[T] 
would only compute [log2(10000)] = 14 lattice ideals. Since the structure of the poset of 
embedded associated lattices can be complicated it is not known if there are better search 
algorithms than random search. 
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